PL-TR-94-2247 


DIELECTRIC  PROPERTIES  OF  MAGNETOPLASMAS 


Francis  S.  Johnson 


University  of  Texas/Dallas 
Box  830688 

Richardson,  TX  75083-0688 


20  November  1989 


Scientific  Report  No.  3 

19950421  052 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


PHILLIPS  LABORATORY 

Directorate  of  Geophysics 

AIR  FORCE  MATERIEL  COMMAND 

HANSCOM  AIR  FORCE  BASE,  MA  01731-3010 


"This  technical  report  has  been  reviewed  and  is  approved  for 
publication" 


This  report  has  been  reviewed  by  the  ESC  Public  Affairs  Office  (PA) 
and  is  releasable  to  the  National  Technical  Information  Service 
(NTIS)  . 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense 
Technical  Information  Center  (DTIC) .  All  others  should  apply  to 
the  National  Technical  Information  Service  (NTIS) , 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  the 
mailing  list,  of  if  the  addressee  is  no  longer  employed  by  your 
organization,  please  notify  PL/TSI,  29  Randolph  Road,  Hanscom  AFB, 
MA  01731-3010.  This  will  assist  us  in  maintaining  a  current 
mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations 
or  notices  on  a  specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATiON  PAGE 


Form  Approved 
OM3  No.  0704-0188 


Pjcl:c  '?30r'"i9  D 


:rden  for  ims  collection  of  information  is  estimated  to  average  i  hour  oer  resportse.  including  the  time  for  reviewing  instruaions,  searching  existing  data  sources. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

20  November  1989 


4.  TITLE  AND  SUBTITLE 

Dielectric  Properties  of  Magnetoplasmas 


3.  REPORT  TYPE  AND  DATES  COVERED 
Scientific  No.  3 


5.  FUNDING  NUMBERS 

PE  61102F 

PR  2310  TA  GL  WU  AA 


3.  AUTMCR(S) 

Francis  S.  Johnson 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
University  of  Texas/Dallas 
Box  830688 

Richardson,  TX  75083-0688 


Contract  F19628-93-K-0008 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


;9.  SFONSORING/MONITORIF^G  AGENCY  NAME(S)  AND  ADDRESS(ES) 
i  Phillips  Laboratory 
i  29  Randolph  Road 
Hanscom  AFB,  MA  01731-3010 

I 

■  Contract  Manager:  David  Anderson/GPIM 


11.  SUr'PLLME.'.VAr.V  NOTES 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


PL-TR-94-2247 


12b.  DISTRIBUTION  CODE 


123.  CISTRiEUTION/AVAlLABILlTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 


13.  ABSTRACT  (Maximum  200  words) 

This  report  discusses  the  dielectric  properties  of  magnetoplasmas  by  first 
considering  a  plasma  at  rest  and  then  causing  acceleration  by  increasing  the 
electric  field,  E,  slowly.  The  polarization  current  is  calculated,  then  gravity 
is  added  which  causes  a  current  to  flow.  It  is  found  that  under  the  J  x  B  force, 
the  plasma  stops  in  time  and  then  reverses  its  direction  of  flow.  Finally,  the 
study  considers  an  alternate  treatment  which  does  not  involve  dielectric 
properties. 


14.  SUBJECT  TERMS 

Dielectric  properties 
Magnetoplasmas 
Electromagnetic  forces 


15.  NUMBER  OF  PAGES 

8 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 


NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  THIS  PAGE  OF  ABSTRACT 

Unclassified  Unclassified  SAR 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std  Z39-18 
298-T02 


DIELECTRIC  PROPERTIES  OF  MAGNETOPLASMAS 
(dielectric,  11/20/89) 


Consider  a  slab  of  plasma  bounded  by  capacitor  plates  and  permeated 
by  a  uniform,  constant,  magnetic  field  in  the  z  direction.  Apply  a  potential 


Oib 


difference  across  the  plates  so  that  E  is  in  the  +y  direction.  Then  the  plasma 
drift  is  in  the  +x  direction  and  Vx  =  E/B.  This  agrees  with  the  first  term  of 
Schmidt’s  equation  (2—166), 

ExB  ,  mEp 


The  second  term  is  in  the  +y  direction  if  E  is  increasing  with  time.  Schmidt 
identifies  this  term  as  the  poiarization  drift  and  notes  that  it  is  associated  with 
the  conversion  of  electric  field  energy  to  plasma  kinetic  energy;  thus  the 

polarization  current  is  in  the  direction  of  E,  and  Jp-E  >  0  for  E  >  0. 

If  we  consider  the  plasma  to  be  initially  at  rest  (and  cold)  and 
accelerate  it  by  increasing  E  slowly,  how  much  electrical  energy  must  be 
supplied  to  the  capacitor  to  do  this?  When  the  electric  field  reaches  the  value 
E,  the  energy  of  the  system  is 

eoE2  ,  Nmv2 

-7i - 1 - ?! -  • 


The  y  component  of  motion  of  the  guiding  center  is  given  by  the  second  term 
in  Schmidt’s  equation,  and  the  displacement  of  the  guiding  center  is 


m  dE  _  mE 
3T  ~  q^  • 


The  energy  extracted  from  the  electic  field  is 


'  -KT-r.  m  dE  NmE2  _  Nmv^ 

dT  “  "T” 

which  is  just  the  kinetic  energy  imparted  to  the  plasma. 


NmE2  Nmv2 


The  polarization  current  is 

m  dE  _  Nm  dE 
p^  ^  cTT  • 
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In  dielectric  terms,  this  contributes  to  the  polarization  P  of  the  medium, 


where 


p  _  NmE  _ 


Nm 

7^ 


£oE. 


Thus  the  electrical  susceptibility  is 


_  Nm  _  Nmv2 

“  7^  -  TTE^ 

and  the  relative  permittivity  is 


fr  —  1  + 


Nmv2 

T^- 


The  electric  field  energy  density,  expressed  in  terms  of  the  dielectric 
properties,  is 


eo^rE^  _  ,  Nmv2 


just  the  total  energy  mentioned  above.  The  kinetic  energy  in  the  medium  is 
recovered  as  electrical  energy  when  the  capacitor  is  discharged. 


It  is  useful  to  look  at  the  magnitude  of  some  parameters  typical  of 
space  plasmas: 


B  =  10  nT  =  10-8  T 
m  =  1.67  X  10-27  kg 
€o  =  8.85  X  10-12 
E  =  Bv  =  10-3  V/m 
P  =  Xe^oE  =  1.68  X  10-7  =  NqAy; 
Proton  gyro  radius  =  mv/qB  =  1 


N  =  107  m-3 

q  =  1.6  X  10-19  c 
V  =  103  m/s 
er  =  1.9  xl07 
Ay  =  1.05  X  103  m 

X  103  in  (smaller  if  thermal  velocity  is  used). 


Thus  the  displacment  due  to  polarization  is  of  the  order  of  100  gyro  radii  or 
more.  The  charge  density  that  must  be  applied  to  the  capacitor  plates  to 
produce  a  velocity  of  103  m/s  is 


CTf  =  fo^rE  =  1.68  X  10-7  C/m2; 

this  would  produce  a  field  of  1.9  x  lo^  V/m  except  for  the  polarization  of  the 
plasma.  This  also  shows  that  eoE2/2  is  generally  negligible  by  comparison 
with  Nmv2/2,  the  importance  of  the  electric  field  in  the  problem  notwith¬ 
standing. 


Remove  the  Capacitor  Plates 

If  a  slab  of  plasma  is  moving  through  a  uniform  B  field  with  velocity  v, 
there  must  be  surface  charges  sufficient  to  produce  an  electric  field  within  the 
plasma  that  corresponds  to  its  drift  velocity.  These  surface  charges  are 
usually  attributed  to  v  x  B  forces.  It  seems  natural  enough  to  call  this 
polarization  of  the  medium,  but  this  inevitably  leads  to  confusion  with  the 
polarization  current  discussed  above,  which  is  in  the  opposite  direction.  To 
avoid  this  confusion,  we  will  refer  to  the  surface  charge  density  as  the  induced 
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surface  charge  density  cri.  whose  magnitude  is  the  same  as  the  free  charge 
density  that  had  to  be  placed  on  the  capacitor  plates  in  the  previous 
discussion.  Note  that  a\  is  positive  on  the  lower  surface  and  negative  on  the 
upper,  whereas  the  surface  charge  density  Cp  associated  with  the  polarization 
is  negative  on  the  lower  surface  and  positive  on  the  upper,  so  the  two  tend  to 
cancel . 

CTi  —  foCrE  =  (1  +  Xe)  foE 

=  Xe  foE 

Thus,  for  space  plasmas,  the  two  cancel  to  within  about  one  part  in  10^,  but  it 
is  just  this  small  difference  between  them  that  is  responsible  for  the  electric 
field  within  the  plasma. 


Now  Add  a  Gravity  Field 

(1)  Start  with  stationary  plasma  between  capacitor  plates,  with  the  plates 
connected  together  so  that  there  is  no  electric  field.  To  this  add  a  gravity  field 


■n 


— >‘3 
06 


in  the  +x  direction;  this  will  produce  a  current  in  the  — y  direction,  with 

J  X  B  =  — Nmg,  or  J  =  -Nmg/B. 


The  current  in  the  plasma  produces  a  force  on  the  plasma  that  cancels  the 
force  of  gravity.  The  current  is  continuous  but  dissipationless,  as  the  potential 
difference  between  the  plates  is  zero. 


(2)  Now,  remove  the  connection  between  the  capacitor  plates.  Positive 
charge  will  accumulated  on  the  lower  plate,  leading  to  the  development  of  an 
E  field  in  the  +y  direction  and  consequent  plasma  drift  in  the  +x  direction. 
There  will  also  be  a  polarization  current  in  the  +y  direction,  and  it  will  almost 
cancel  the  current  flow  caused  by  the  presence  of  gravity.  We  can  consider  the 
current  flow  due  to  gravity  as  ^1  accumulating  on  the  capacitor  plates,  with 
the  force  of  gravity  being  totally  cancelled  by  the  current.  But  then  we  must 
recognize  that  there  is  a  polarization  current  in  the  plasma  that  accelerates  the 
plasma  in  the  +g  direction,  and  that  this  acceleration  is  smaller  than  g  by 
only  about  1  part  in  10^. 


What  must  actually  happen  is  that  the  two  currents  largely  cancel  one 
another,  and  the  plasma  is  nearly  in  free  fall.  A  very  small  part  of  the  work 
done  by  the  gravity  field  produces  electric  field  energy  eoE2/2,  and  the 
predominant  remainder  goes  into  the  production  of  kinetic  energy. 

(3)  It  is  not  necessary  to  imagine  the  presence  of  the  capactitor  plates  at 
^1.  A  slab  of  plasma  in  a  uniform  field  B  will  accelerate  at  almost  the  rate  g, 
while  accumulating  net  charge  density  at  its  surface  a  =  a\  —  a^,  electric  field 
in  the  +y  direction  E  =  c/eo,  and  electric  field  energy  density  eoE2/2  in 
addition  to  its  kinetic  energy  density. 
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For  typical  space  plasmas,  we  noted  that  charges  were  displaced 
something  like  100  km  due  to  polarization.  In  the  case  of  a  plasma  slab 
accelerated  by  gravity  to  100  km/s,  the  net  surface  charge  density  is 
approximately  cTp/^e,  and  the  charge  displacements  required  to  produce  the 
electric  field  is  only  about  1  cm. 

(4)  Now  let  g  act  to  slow  down  a  slab  of  plasma  that  is  initially  moving 
with  velocity  Vq  in  the  +  x  direction.  From  the  preceding,  the  net  surface 

n - 


o  ^  ^ _ _ _ 

charge  initially  is  cr  =  CqE  =  eoBvo.  This  can  be  thought  of  in  terms  of 
(J  =  (Ti  -  o-n,  or  (7i  =  <7  +  0-p,  whoro  (7p  =  P  =  XeCoE;  i.e.,  ai  =  eoE  +  Xe^oE 
=  eoCrE.  (This  is  just  the  cnarge  that  could  be  withdrawn  from  the  capacitor 
plates  in  an  earlier  discussion.) 


Now  the  g  field  causes  a  current  to  flow  in  the  +y  direction,  and  the 
J  *  B  force  is  in  the  +x  direction,  again  opposing  gravity.  Suppose  for  the 
sake  of  argument  that  it  exactly  cancels  gravity.  The  current  flow  now  acts  to 
decrease  the  surface  charge,  and  would  reduce  it  to  zero  in  a  time 

rr-ll  —  ^  ^ 

“  Ning^ 

(at  which  time  the  E  field  would  be  reduced  to  zero  and  the  plasma  would  be 
stationary,  but  we  must  still  show  this). 


The  tendency  of  the  current  Nmg/B  to  decrease  the  surface  charge  is 
largely  counteracted  by  the  decreasing  polarization  of  the  medium.  As  the 
charge  eo^rE  is  delivered  to  the  upper  boundary  by  the  current  that  cancels 
the  force  of  gravity  on  the  plasma,  the  polarization  charge  on  the  upper 
surface  is  reduced  from  Xe^oE  to  zero,  and  the  depolarization  current 
associated  with  this  accelerates  the  plasma  in  the  — x  direction,  producing  an 


acceleration  almost  equal  to  g.  It  is  where  Xe  -  10^-  The  fact  that  the 

plasma  does  not  slow  down  quite  as  fast  as  it  would  under  action  of  gravity 
alone  reflects  the  fact  that  dectric  field  energy  eoE2/2  is  eliminated  in  the 
process,  contributing  to  the  maintenance  of  velocity  in  the  +x  direction. 


Under  action  of  the  accderation  produced  by  the  Jp  x  B  force,  the 
f  f  FR 

plasma  stops  in  time  ° —  ,  and  then  reverses  its  direction  of  flow  and  the 

direction  of  the  electric  field  within  it. 
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Alternate  Treatment  Not  Involving  Dielectric  Properties 

It  is  probably  preferable  to  analyze  the  response  of  a  slab  of  plasma  to  a 
gravitational  field  without  invoking,  or  developing,  the  dielectric  properties  of 
the  magnetoplasma.  As  before,  consider  a  slab  of  plasma  in  the  xz  plane, 

I 


I  ^ 

I _  -I 

Gg-  ^ 


with  gravity  acting  in  the  +x  direction.  Consider  a  current  J  in  the  +y 
direction;  J  will  turn  out  to  be  negative  so  as  to  produce  a  J  x  B  force  on  the 
plasma  that  opposes  gravity.  Later  on,  we  will  solve  for  J.  Newton’s  Law 
gives  for  the  acceleration  of  the  plasma 

Nma  =  Nmg  +  JB,  or  a  =  g  +  JB/Nm. 

The  drift  velocity  must  be  consistent  with  the  electric  field  that  results  from 
the  surface  charges  produced  by  current  J.  The  surface  charge  on  the  upper 
surface  of  the  slab  is 


(J  =  /  J  dt 

and  E  =  -a/co  (fr  =  Ij  as  no  dielectric  properties  are  involved  in  this 
analysis).  For  consistency,  it  is  necessary  that  v  =  E/B  =  /  a  dt,  or  that 


Therefore 

and 


-J/eoB  =  g  +  JB/Nm  . 

T If  ^  _  Nmg  1 

“"^NmvVeoB2  +  1  ’ 

^  =  (Nmg/eoB2) 

=  6  M„„,i  1 


+  1 


1  +  CoB2/Nmv2  ■ 

This  is  smaller  than  g  by  about  1  part  in  10^.  It  is  the  same  expression  as 
g  derived  earlier  making  use  of  the  dielectric  properties  of  the  plasma. 
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